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Abstract. Using the method presented in Erspamer & North ( 2002, Paper I hereafter), detailed abundances of 140 
stars are presented. The uncertainties characteristic of this method are presented and discussed. In particular, 
we show that for a S/N ratio higher than 200, the method is applicable to stars with a rotational velocity as 
high as 200 km s" 1 . There is no correlation between abundances and V svai, except a spurious one for Sr, Sc and 
Na which we explain b y the small number of lines of these elements combined with a locally biased continuum. 
Metallic giants (Hauck, 1986) show larger abundances than normal giants for at least 8 elements: Al, Ca, Ti, Cr, 
Mn, Fe, Ni a nd Ba. The anticorrelation for Na, Mg, Si, Ca, Fe and Ni with V smi sugge sted by Varenne and 
Monier ( 1999] ) is not confirmed. The predictions of the Montreal models (e.g. Richard et al., 2001) are not fulfilled 
in general. However a correlation between and logp is found for stars of 1.8 to 2.0 Mq. Various possible 

causes are discussed, but the physical reality of this correlation seems inescapable. 
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1. Introduction 

A-F type stars present a large number of chemical pecu- 
liaritie s: Ap, Am, A Bootis, p Puppis (see Gray & Garrison 
|l989b , for a discussion of the 8 Del stars and a definition 
of the p Pup stars). Radiative diffusion is generally in- 
voked most of these peculiarities, and the theory of this 
mechanism, developed mostly by Michaud and collabora- 
tors (Michaud 1991) seems indeed largely successful, even 
though some adjustments seem to remain necessary. For 
instance, stellar wind was added to radiative diffusion in 
order to make more realistic predictions of abundances of 



Am stars (Michaud et al. |1983[ ). 

A little explored category of chemically peculiar stars 
in this range of spectral types is that of metallic giants, 



discovered by Hauck ( |l986[ ) on the basis of t cxtscGcneva 
photometry and further studied by Berthet ( 1990 , 1991 ) 
using spectroscopy. They are A-F stars leaving the main 
sequence and having the Geneva photometric parameter 
Am.2 > 0.013 (see Hauck, 1973 for the definition of Ama). 
Berthet (199C, 1991) confirmed that these stars have metal 
abundances higher than solar ones, as suggested by their 
large Am2 values. However, he analyzed only metallic gi- 
ants and relied on published abundances of normal giants 
for comparison. Berthet found that the abundance pattern 
of metallic giants is similar to that of Am stars, except 
for calcium and scandium, which are notoriously under- 
abundant in the latter, while they are roughly normal or 
slightly overabundant in the former. On this basis, Berthet 



(1992) proposed that metallic giants might be evolved Am 
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stars, an idea theoretically justified by the following sce- 
nario: in the beginning of the life of an Am star, Ca and 
Sc sink and accumulate below the line forming region be- 
cause of their too small radiative acceleration. During stel- 
lar evolution on the main sequence, however, the effective 
temperature decreases and the outer convection zone de- 
velops and thickens, until it reaches the depth where Ca 
and Sc have accumulated and dredges them up, so that 
their surface abundances becomes normal again. Kiinzli 
& North (1998) had challenged this scenario, arguing, in 



2 



D. Erspamer and P. North: Abundances of 140 A-F stars with ELODIE and CORALIE 



particular, that the V sini distribution of metallic giants 
is not compatible with that of Am stars, especially evolved 
ones. The fact that metallic giants are more evolved and 
more massive than normal ones (North et al. 1998) does 
not fit Berthet's scenario either. 

The initial purpose of this work was to focus our at- 
tention on metallic giants by comparing their abundances 
with those of normal giants, in order to better define their 
characteristics and determine how far they form a well 
defined spectroscopic group. However, since Berthet's sce- 
nario remains unsatisfactory, the question of their origin 
remains open, so we decided to examine these stars in a 
more general context and try to understand what are their 
progenitors. So we judged important to study not only gi- 
ant A-F stars (normal and metallic) but also their possible 
progenitors, i.e. normal A-F type stars of similar masses 
on the main sequence. To this purpose, we chose to ob- 
serve a sample of A-F stars that cover the HR diagram as 
homogeneously as possible, in order to map the abundance 
pattern in the whole region of interest. In addition to the 
homogeneous coverage of the HR diagram, it is important 
to ensure also the homogeneity of the method of spec- 
tral analysis. In the latter sense, homogeneity is crucial, 
since it is very delicate to compare abundance analyses 
done by different authors. There is always at least one as- 
sumption that differs. The most common problem is the 
choice of atmospheric parameters, T e ff and logg. Indeed, 
different techniques, using different photometric systems 
or spectroscopy, are used to derive atmospheric parame- 
ters. Abundance differences arising from these techniques 
are generally not very important since they remain within 
error boxes, but they significantly increase the scatter in 
the results. Another common difference is the choice of at- 
mosphere model theor y. In t he field of A-F stars, ATLAS9 
models from Kurucz ( 1993 ) are the rule. However, these 
models exist with and without overshooting, and both are 
used depending on the authors. Finally, methods to de- 
rive the abundances differ. Some people use equivalent 
widths while others prefer synthetic spectra adjustments. 
The oscillator strengths and other line parameters may 
come from different sources and these sources may vary 
from one element to the next. The adopted solar iron 
abundance may also differ from one author to the other. 
Therefore, even for iron, a large scatter is found. For ex- 
ample, even for a very bright star such as Procyon, [Fe/H] 
values vary from -0.18 to 0.05 (Cayrel de Strobel et al. 
2001 ). All this shows how well justified is the extensive 
and homogeneous analysis of a large sample of A-F type 
stars that we present in this paper. 

In Sect. |, we describe the observations and the re- 
duction procedure, while Sect. [| gives more details on 
the definition of the sample. The determination of atmo- 
spheric parameters is described in Sect. |J, and sources of 
uncertainties on the derived abundances are discussed in 
Sect. ||. The results are presented in Sect. ||and discussed 
in Sect. @. 



2. Observations and data reduction 

The spectra used in this work were obtained with two 
almost identical spect rograp hs, ELODIE and CORALIE 
(see Baranne et al., 1996, for technical details about 
ELODIE and Queloz et al., |200l[ for CORALIE). Both 
spectrographs have been built as a joint project between 
the Geneva Observatory and the Observatoire de Haute- 
Provence. ELODIE is attached to the 1.93 m telescope of 
the Observatoire de Haute-Provence whereas CORALIE 
is attached to the Swiss 1.2 m Leonard Euler telescope at 
the ESO la Silla Observatory. The optical design of both 
spectrographs is the same. The major difference is the 
better resolving power of CORALIE (50 000 compared to 
42 000 for ELODIE) that comes from a slightly different 
optical combination at the entrance of the spectrograph 
and the use of a 2k x 2k CCD camera (against a lkx lk for 
ELO DIE) with smaller pixels (15 /im vs 24 /mi) (Queloz 
et al. |1999| ). 

Two observing runs were held in June and August 1999 
in Haute-Provence, and one in July 2000 in Chile for a 
total of 23 nights (10 and 13 respectively). Target stars 
were bright ones, because we wanted a high signal to noise 
ratio (S/N > 200). The total number of observed stars is 
140, after elimination of all spectroscopic binaries or stars 
presenting strange cross-correlation functions due e.g. to 
pulsations. 

Data reduction was done as presented in Paper I. 
As the optical design of CORALIE is the same as the 
one of ELODIE, only a few changes in the parameters 
used to characterize the CCD were needed to adapt the 
IRAF reduction procedure which was initially designed 
for ELODIE spectra. However, after order extraction, the 
different orders merge perfectly without scaling, contrary 
to the case of ELODIE spectra. The whole spectrum is 
obtained by a simple concatenation of consecutive orders 
with a detection of the overlapping region. As the S/N 
ratio is better in the red part of the order than in the blue 
one, the flux of the first order is retained for 3/4 of the 
overlapping region, and the flux of the following for the 
remaining 1/4. The subsequent operations are done using 
the same procedure as described in Paper I. 



3. Sample 

The mass range is 1.3 to about 2.8 Mq, encompassing that 
of most metallic giants (1.7 to 3 M Q ). As we wanted to ob- 
serve as many stars as possible during the observing runs 
at our disposal, our choice naturally went to the brighter 
stars. Two additional requirements were used to define the 



sample: availability of Geneva photometry (Golay 1980; 
Rufener & Nicole t |1988| ; Cramer |1999| ) and of Hipparcos 
parallaxes (ESA 1997 ), in order to be able to determine 



the atmospheric parameters with the method presented 
in Sect. || Known binaries were exluded too. Finally, we 
chose to keep only stars with V^sinz < 150 km s _1 or 
that have no determination in the bright star catalogue 
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(Hoffleit & Jaschek 1991 ). These were the only restrictions 
before observing. 

During the observing runs, as our main targets were 
metallic giants, we focused a little more on them and on 
their progenitors, i.e. on stars falling in the mass range 
1.7 to 2.5 M Q . Low mass stars had been included in the 
sample essentially because the early Montreal models (see 
Tur cotte et al. ( |1998| ), Richer et al. fl200C| ) and Richard et 
al. (2001) and reference therein) did not consider masses 
larger than 1.5 M Q . 

As observations were done with spectrographs de- 
signed to get radial velocities, the cross-correlation peak 
was routinely available, we used it to exclude stars pre- 
senting a "dubious" peak, i.e. an asymmetric one that 
may come either from binarity or from pulsation, or a very 
broad one resulting from Tasini larger than 200 km s . 

The final sample is composed of 140 stars, out of which 
35 are giants (18 metallic and 17 normal), 10 ar e Am and 
6 are Ap. Spectral types are taken from Hau ck (1986) for 
stars i n common, and from Gray & Garrison ( 1987 , 1989a , 



log- 



L 



-0A(M y - 4.75 + B.C. 



1989b ) otherwise. The chemically peculiar Am and Ap 



stars are included only because their positions in the HR 
diagram correspond to the studied region. It is interesting 
that without taking care of this detail, the proportion of 
Am and Ap stars, relatively to stars of the main sequence, 
is representative. 



4. Atmospheric parameters determination 

Atmospheric parameters were determined from Geneva 
phot ometry for T cS and from HiPPARCOS d ata (E SA, 
1997| ) for log 5. The calib code (Kiinzli et al., jl997| ) was 
used to derive effective temperature from textscGeneva 
photometry. Although this code also provides log (7, we 
prefer to determine that value from HiPPARCOS data. 
All stars of the sample are bright and have an accurate 
trigonometric parallax measurement by HiPPARCOS. They 
all have a parallax larger than 7.4 mas (i.e. are closer than 
135 pc except two, HD 175510 and HD 187764 that are at 
166 and 171 pc respectively). Therefore reddening was ne- 
glected [| Bolometric luminosity was calculated using par- 
allax, bolometric correction interpolated in Flower's table 
(|1996|) and M bo i = 4.75 for the Sun. The bolometric lu- 
minosity is then given by 



1 Although this is completely justified for the 121 objects 
closer than 70 pc on the basis of uvby/3 photometry, a closer 
examination requested by the referee shows that the remaining 
objects are slightly more reddened on average. In particular, 4 
stars have E(b - y) > 0.020: HD 115604 (E(b - y) = 0.025), 
HD 181333 (0.035), HD 200723 (0.023) and HD 214441 (0.027), 
implying that T e g might be underestimated by ~ 300 K in 
the worst case. Among the 121 nearby stars, 8 have also 
E(b — y) > 0.020, but they constitute the queue of a fairly 
gaussian distribution, 6 of them being closer than 50 pc. The 
colour exc esses were es timated using the usual calibrations by 
Crawford ( 1975 , 1979 ), incorporated in a code kindly made 
available by Dr. C. Jordi. 



The logg value was then interpolated in the grids of 



models of Schaller et al (1992), Schaerer et al. (1993a), 
Charbonnel et al. (1993) and Schaerer et al. ( 1993b[) ac - 
cording to the method described by North et al. (1997). 
Finally, an ATLAS9 atmosphere model was computed 
without overshooting (Castelli et al., |199€ ) with tempera- 
ture rounded to 25K and log g to 0.05 dex using the version 
of ATLAS9 adapted by M. Lemke0 for UNIX system. 

Derived parameters of our sample stars are presented 
in Table 5 available at CDS. It contains HD, T e g, logg, 
Vrad, Usini and £ m i cr o- The last three values are those 
derived with our method. Fig. |l| shows their position in 
the HR diagram. Stars cover the whole main sequence 
and extend slightly beyond. 




Fig. 1. HR diagram of all studied stars with some evo- 
lutionary tracks (Schaller et al, 1992). ■ : stars with 
Am 2 < 0.013; A : stars with Am 2 > 0.013 



5. Uncertainty estimation 

5.1. Errors resulting from wrong atmospheric 
parameters 

In order to check the influence of erroneous atmospheric 
parameters (T e g ±150-ftf , logg ±0.2 dex) on derived abun- 
dances, simulations were done on a synthetic spectrum. 
These values do not represent the internal error of the 
methods used to derive these parameters, but rather the 



2 http:/ /www.sternwarte. uni-erlangen.de/ftp/michael/ 
atlas-lemke.tgz 
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Table 1. Difference between adjusted and reference abun- 
dances in hundredth of dex. The last column is the result 
of the adjustment done with the correct atmosphere model 
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-0 


Mg 
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-3 
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-0 


Al 
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-1 





-1 


Si 
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-4 





-1 


-1 
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-2 
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-4 


-1 


Ca 
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-8 


-1 


1 


-0 


Sc 
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-5 
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-8 


-1 


Ti 


7 


-7 


4 


-4 


-0 
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-7 
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-4 


-0 


Cr 
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-6 
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-2 


-0 


Mn 


8 


-8 











Fe 


8 


-8 


-0 


-0 


-0 


Co 
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-8 


1 


-1 





Ni 


7 


-9 


-1 


-0 


-1 


Cu 
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-10 


-0 


-0 


-0 


Zn 


8 


-8 











Sr 


10 


-9 


2 


-1 


-0 


Y 


6 


-6 


7 


-6 


-0 


Zr 


5 


-5 


7 


-7 


-0 


Ba 


9 


-9 
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-3 


-1 


La 


8 


-5 


8 


-5 


1 


Ce 


7 


-7 


6 


-6 





Nd 


9 


-8 


7 


-6 


-0 


Sm 


8 


-7 


6 


-6 


-0 


Eu 


8 


-8 


6 


-6 






scatter in the estimations done with different methods 
(photometric or spectroscopic). The synthetic spectrum 
was computed with a model atmosphere of T c ff = 6950 K, 
log<7 = 4.0 and solar chemical composition. This model 
corresponds to a typical star of our sample lying in the 
middle of both temperature and gravity ranges. The other 
parameters, namely radial, microturbulent and especially 
rotational velocities, are V ra a — km s" 1 , £micro = 1.5 
km s _1 and V sini = 10 km s . The latter was chosen 
so that the code returns the input abundances within 
±0.01 dex when abundances are estimated with the input 
atmosphere model, the S/N ratio being infinite except for 
rounding errors. Thus all variations in derived abundances 
can be attributed to errors on T e ff or logg, even though 
^ad, £micro an d V sin i were adjusted, as well as the abun- 
dances. The input abundances are solar. The results are 
presented in Fig. || and listed in Table |l| 

An erroneous T c g leads to an error in abundances < 
0.1 dex with almost all elements affected in the same way 
except C, O, and S which react in the opposite direction 
(see Fig. || b and c) . 



It is interesting to note that an erroneous log g has al- 
most no effect on abundances of iron peak elements (Cr to 
Zn), and leads to an error < 0.1 dex for the other elements. 

When combining both errors, the total error is always 
< 0.15 dex, as elements that are very sensitive to temper- 
ature are not very sensitive to gravity and conversely, at 
least for the adopted temperature and gravity. 

5.2. Errors resulting from wrong continuum 
adjustment 

Another source of error comes from bad continuum ad- 
justment. When rotational velocity increases, it becomes 
harder to find continuum points. The use of echelle spec- 
tra is helpful because the whole spectrum from 3900A to 
6820A is at our disposal. The normalization method is 
presented in Paper I but can be summarized as follows : 
the spectrum is divided into 6 parts, with one common or- 
der between parts, and each part is normalized separately. 
Each part has a length of more that 400A and was chosen 
so that there are some continuum points near the borders. 
The continuum is fitted with the continuum function of 
IRAF whose parameters were defined using slow rotators, 
mainly the Sun and Vega. 

Even when working with high S/N ratios, the main 
problem for continuum estimation is rotational velocity. It 
is relatively easy to normalize a spectrum for a star with a 
V smi <40 km s — 1 , because lines are not too broad. When 
rotational velocity increases, lines become broader and it 
is hard to find lines that are not blended. Blends become 
so important that continuum points vanish. Selecting wide 
parts allows to preserve continuum points even for large 
rotational velocities. 

Simulations were done on synthetic spectra having the 
same parameters as in previous section, but with V sini 
= 10, 50, 100, 150 and 200 km s _1 . Moreover, in order to 
simulate as closly as possible real spectra, gaussian noise 
was added to obtain a reference spectrum with S/N = 
200, which corresponds to the minimum observed value. 
Finally, the spectra were scaled by a factor 0.99 and 1.01, 
so that the continuum level is respectively too high (by 
a factor 1.01) and too low (by a factor 0.99). For larger 
differences, it becomes obvious that the continuum is false. 
Results of these test are presented in Fig. || and in Table |[ 

It is interesting to note that for the correct continuum 
definition, almost all adjusted abundances are well within 
0.1 dex of their input values, even for high V sin i (see cen- 
tral column of Fig. [|). Moreover the adjustment with the 
correct value allows to select "sensitive" elements, such as 
Al, Si, S, V, Cu, Y, and Nd. These elements have very few 
and very weak lines, generally blended with much stronger 
lines. 

Another interesting point is the shape of the abun- 
dance pattern obtained with a scaled continuum. This 
shape remains very constant when the rotational velocity 
increases, but its amplitude increases with this velocity. 
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(a) T eff = 6950 ; log g = 4.0 



(b) T eff = 6800 ; logg = 4.0 



(c) T eff = 7100 ; lo gS r = 4.0 



C 14a Mg A] Si S CaScTi V CrMnFeCo Ni CuZnSr Y Zr Ba La CeNdSmEu 



C NaMg Al Si S CaSc Ti V CrMn Va Co Ni Cu'/ri Sr Y Zr Ba La CeNdSmEu 
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(e) Teff = 6950 ; log 5 = 3.: 



(f) Teff = 6950 ; logg = 4.2 



Fig. 2. Abundances with respect to the Sun for different values of T e g and \ogg. (a) Abundances obtained with the 
correct atmosphere model, (b) to (e) Abundances obtained with the specified parameters. Error bars represent the 
sigma of the average abundance value of the 7 parts of the spectrum (see Paper I) . 



Finally, this simulation can be interpreted as yielding 
the maximum possible error, because all adjusted parts 
were scaled by one and the same factor. In reality, even for 
V sini = 150 km s _1 , only the first two parts (A < 4900A) 
can be affected by the lack of continuum points, the 5 re- 
maining parts being much easier to normalize. Therefore 
it is possible to conclude that the influence of the normal- 
ization on real spectra will be well below the values listed 
in Table |. 

6. Results 

All results are listed in Tables 5 and 6 available at the 
CDS. The first table lists, for each star, the radial ve- 
locity, the microturbulence and the projected rotational 
velocity. The second one lists the abundances. In 12 cases, 
there are two spectra per star - hence two lines per star 
in Tables 5 and 6 - either because the same object was 
observed once with ELODIE and once with CORALIE (6 
cases) , or because the star was observed twice with a given 
instrument (6 more cases). 

6.1. Abundances 

As the sample covers a large temperature domain, it is 
difficult to summarize the results. For early A-type stars, 
only a few elements show lines in the spectrum while for 



F-type stars, it is possible to derive abundances for 28 
elements. Table 6, available at CDS, contains the detailed 
abundances for each star and the number of lines used for 
the determination. 

24 stars of the sample are common with the cata- 
logue of [jf] done by Cayrel de Strobel et al. (2001). 



Comparisons for T e g, logg and [^1 are illustrated in 
Fig. |j. The agreement is very good. It is intersting to no- 
tice that 12 stars come from a single, homogeneous source, 
Balachandran (1990). This sample alone shows a slight 
shift in Toff of less than 100 K and log g values which are in 
good agreement with ours, [^jf] values are in good agree- 
ment, except for two cases where the difference amounts 
to about 0.3 dex. 



6.2. Microturbulence 

Microturbulence is an ad hoc parameter which relates to 
turbulent motions on scales smaller than the mean free 
path of a photon. As commonly done in practice, we as- 
sumed it is constant with optical depth even though it 
should be allowed to vary in principle. Our estimates are 
based on the the four bluer spectral segments only, since 
there are too few appropriate lines in the three redder seg- 
ments to constrain the adjustment; the microturbulence 
value is an average of the four individual fitted values. 
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Element Element Element 



(a) Vsini = 10 ; S/N = 200 ; (b) V sini = 10; S/N = 200 ; cont = 1 (c) Vsini = 10 ; S/N = 200 ; 
cont = 1.01 cont = 0.99 



C C NaMg A. Si S CuScTi V CrMnFe Co Xi CuZr. Sr Y ZrDaLaCeXcSmEu 



C NaMg Al Si S Ca Sc Ti V CrMnFeCo Ni CuZnSr Y ZrBaLa CeNdSmEu C NaMg Si S CaSc Ti V CrMnFe Co Ni CuZn Sr Y Zr Ba La CeNdSmEu 



TO 




(d) Vsini = 50 ; S/N = 200 ; (e) V sin i = 50; S/N = 200 ; cont = 1 (f) Vsini = 50 ; S/N = 200 ; 
cont = 1.01 cont = 0.99 



C NaMg A". Si S Ca Sc Ti V CrMnFe Co Xi CuZr. Sr Y Zr 3a Li. CeXcSirEu ; CO XuMcM Si S Ca Sc Ti V CrMnFe Co Xi Ci.iZr. Sr Y ZrBaLa CeXcSmEu 



C NaMg Si S CaSc Ti V CrMnFe Co Ni CuZn Sr Y Zr Ba La Ce NdSmEu 




(g) Vsini = 100 ; S/N = 200 ; (h) Vsini = 100 ; S/N = 200 ; (i) Vsini = 100 ; S/N = 200 ; 
cont = 1.01 cont = 1 cont = 0.99 



:: N.iMs A :-i S Cii Sc Ti V CrMnl oCo \i Cu'/r. ;-v V Zr ia I ,a Co XtlSnil ',u V \a\'j;A Si S C'aS; li Ci iin i :i Co V Cii /it Sr V /■ in I ,;: Cc \i:S;i: in 




C NaMg Si SCaScTi V CrMnFe Co Ni CuZn Sr Y Zr Ba La Ce NdSmEu 



V 



(j) Vsini = 150 ; S/N = 200 ; (k) Vsini = 150 ; S/N = 200 ; (1) Vsini = 150 ; S/N = 200 ; 
cont = 1.01 cont = 1 cont = 0.99 



Fig. 3. Abundances with respect to the Sun for different continuum heights. Error bars represent the sigma of the 
average abundance value of the 7 parts of the spectrum (see Paper I) . 
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Table 2. Difference between adjusted and reference abundances in hundredth of dex for different rotational velocities 
and different continuum height. 
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Because of the large size and homogeneity of our sample, 
the relation between microturbulence, effective tempera- 
ture and luminosity is worth examining. 

The values are reported in Table 5 available at CDS. 
Fig. U shows the relation between T c g and £ m i cr o for reli- 
able £ m icro estimates, defined as those for which the stan- 
dard deviation of the 4 individual fitted values is smaller 
than 0.5 km s . This relation is compared with the av- 
erage, hand-drawn one presented in Coupry & Burkhart 
( |1992p , which is the continuous curve in Fig. §. The agree- 
ment is good for cool stars but, beyond T e g = 7400 K, the 
relation of Coupry & Burkhart becomes an upper limit to 
our values. 



The relation of Nissen ( |1981| ) , which is linear in both 
T e ff and logy, is represented by the three straight lines in 
its domain of validity. The agreement seems not too bad 
for low luminosity stars, but our data seem to indicate 
that all stars share the same £ m i cro at low temperature, 
whatever their log <?, while evolved stars separate from the 
unevolved ones very rapidly as T e g increases. The T e ff de- 
pendence of £ m icro is much steeper in our case than in 
Nissen's, and the slope seems to increase with decreas- 



ing logg. Qualitatively, however, we confirm that evolved 
stars have higher £ m i C ro in this range of T e g. This is also 
in qualitative agreement with the results of Gray et al. 
(2001) for stars cooler than A8. 



Strangely enough, the trend reverses for T e ff between 
8000 and 8800 K, evolved stars having significantly lower 
Cmicro on average than unevolved ones. The latter trend is 
largely due to the group of stars with T c s larger than 8000 
K and £ m i C ro larger than 2.4, which correspond to stars 
with mass of 1.8 to 2.0 M Q that present iron and nickel 
overabundances (see Sect. 7.3 for details). These stars are 
very close to the Zero Age Main Sequence (ZAMS) in the 
HR diagram. This behaviour was already noted by Gray 
et al. ( |2001 ). Whether £ m i cro depends on [M/H] remains 
an open question, and more high quality data would be 
needed to better understand what is happening in this 
diagram. Let us notice that the observed overabundance 
of iron cannot be the result of a bias in the estimate of 
Cmicro, since only an underestimate of it would lead to 
an overestimate of iron abundance; what we see is, on 
the contrary, rather high £ n 
stars. 



values for these metal rich 
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Fig. 4. Comparison between T e g, \ogg, [4^1 from this work and from the [4^1 catalogue of Cayrel de Strobel et al. 
(2001). Filled circles distinguish Balachandran's (1990) contribution, open circles representing determinations by other 
authors. The straight line is not an adjustment but represents the equality of both values 
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Fig. 5. Relation between T e g and £ m icro- A for stars with 
a reliable £ m i C ro determination (see text): stars with log g 
> 4.25. ■ : stars with 4.0 <log g < 4.25. □ : stars with 
3.75 <logg < 4.0. x : stars with 3.5 <log# < 3.75. O : 
stars with log g < 3.5. The continuous line is the relation of 
Coupry & Burkhart (1992), while the thick straight lines 
are Nissen's (1981) relation for the three indicated logy 
values. 



It is difficult to draw any conclusion about a possible 
relation between £ m i cro and log g for the stars hotter than 
8800 K because their number is too small. In addition, 
the £ m icro estimates are less reliable than for cool stars 
because of the much smaller number of lines. 



7. Discussion 

7.1. Rotational velocities 

7.1.1. Comparison with other works 

Fig. ^ illustrates the comparison of our determinations 
of Vsin i with works of Uesugi & Fuku da (1982), Abt & 



Morell ( |1995D , and Royer et al. (|2002[) . Our values are 
slightly larger than those of Abt & Morell, especially for 
V sin i larger than 100 km s^ 1 . However, the trend is in the 
other direction when comparing with Uesugi & Fukuda. 
Finally, comparison with Royer et al. is excellent if we 
neglect one point (^150 km s" 1 for our estimate). This 
point corresponds to HD 12311. This star is still clas- 
sified as a dwarf in Simbad even though it has the at- 
mospheric parameters of a giant. This was already noted 
by Hauck ( 1986) u sing Geneva photometry, and Gray 
& Garrison (1989b) give a spectral type of F0 Ill-IVn. 
Royer et al. ( 2002 ) use Fourier transforms of several line 
profiles to estimate V sini. This technique is very efficient 
for stars with well separated lines. In cooler stars, it be- 
comes harder to find isolated lines to analyze. HD 12311 
is precisely one of the coolest stars of Royer et al. (2002). 
In this case, all lines of their list except one were unus- 
able. Therefore the difference can be explained by a poor 
precision of their determination for this star. In this com- 
parison, there are less points because Royer et al. focus on 
A type stars. 

We conclude that our determination of V sin i are prob- 
ably even better than those of Royer et al., because we 
make use of all lines present in the spectra, whether they 
are blended or not. On the other hand, we neglect macro- 
turbulence, so our V sin i values are slightly biased towards 
large values in very slow rotators (V sin i < 8 — 10 km s _1 ). 
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Fig. 6. Comparison between Vsini values from this work and Uesugi & Fukuda ( |1982| ), Abt & Morcll fll995| ) and 
Royer et al. (2002). The straight line is not an adjustment but represents the equality of both values. 



7.1.2. Systematic abundance error due to rotation 

When dealing with large rotational velocities, one can eas- 
ily underestimate the continuum level. That leads to un- 
derestimated abundances because lines seem shallower. 

The apparent abundances of three elements are 
strongly correlated with V sin i. These elements are Sr, Sc, 
and Na. This is illustrated in Fig. |7|, where iron was also 
included to show the behavior of an uncorrelated element. 
The correlations involving Sr, Sc and Na are spurious, 
however, for the following reasons : 

Strontium: This element shows the most significant cor- 
relation (correlation coefficient r = —0.61). The problem 
arises when only one line, Srii A4215.519, is visible. This 
happens mainly in the hotter stars of the sample. The 
fact is that the continuum is slightly underestimated in 
this region. When rotational velocity increases, this small 
error leads to large underestimates of the abundances. As 
there is no other line of Sr for those stars, this error can- 
not be corrected. In Fig. |7[ almost all stars with Vsini 
> 150 km s _1 are of spectral type A, which reflects the 
well-known drop of rotational velocity for spectral types 
F and later. It is for these relatively hot stars that only 
one Srii line is used, and in this case, the result is much 
poorer than in the simulation illustrated by Fig. which 
involves an F star where two Sr lines are used. 

This is the only negative correlation for all adjusted 
elements. It is very interesting to illustrate the accuracy 
of the continuum adjustment, and the stability of this 
method. Indeed, an underestimate of the continuum leads 
to a negative correlation. The fact that only one element 
shows such a behaviour is a good justification for our 
choice of continuum adjustment. 

Scandium: This element yields the largest positive corre- 
lation (r = 0.5). It should be remembered that during 



the analysis, spectra are cut into 7 parts, each part is ad- 
justed, and finally the different estimates are put together 
(see Paper I for details). In the case of Sc, there is a large 
scatter between abundances of different spectral ranges. 
Moreover, since Sc is ionized, the adjustment is very sen- 
sitive to gravity (see Sect. 5.1). However, there is no corre- 
lation between logg and V sini. Unfortunately, the reason 
of this behaviour was not found, but it is probable that 
there are several causes adding up. 



Sodium: The last significant correlation (r = 0.43) is 
found for Na. For this element, the reason was found to be 
the telluric lines. The Na I doublet (A 5889 and A 5895) 
lies within telluric lines. Even if the adjustment is not 
very sensitive to these lines, a small error can lead to a 
large abundance overestimate, especially for fast rotators, 
because these lines are broad and saturated. 

There is no significant correlation between other ele- 
ments and Vsini. 



7.1.3. Correlation between rotation and true 
abundances 

The decrease of |"-^ ] wit h increasing Vsini found by 
Takeda & Sadakane (1997) and confirmed by Varenne & 



Monier (199E) is firmly denied by our results. Fig. || shows 
clearly that there is no dependence of iron or nickel abun- 
dance with Vsini. In our case, there are 60 stars with 
Vsini > 100 km s^ 1 while there was only 10 in Varenne 
& Monier and 5 in Takeda and Sadakane. It is important 
to notice that the anticorrelation was found for stars of the 
Hyades, but, a priori, there is no reason why only these 
stars should be affected. The result of Takeda et al. (1997) 



is discussed in more details in Sect. 7.3 
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Fig. 8. Abundances of Fe and Ni as a function of V smi for fast rotators. The dashed line represent the anticorrelation 
proposed by Varenne & Monier ( 1999| ). 



Table 3. Characteristics of the observed metallic giants. Table 4. Characteristics of the observed normal giants. 
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7.2. Metallic giants 

Fig. U shows the position of all giant stars of our sample 
in the HR diagram. There are 35 giants, out of which 18 
are metallic. Their characteristics are reported in Tables || 
and|. 

Abundances of observed giants are illustrated in 
There are no very clear-cut peculiarities such 
underabundance of Ca and/or Sc in Am stars. 
However, abundances of some elements (Al, Ca, Ti, Fe, 
Ba for example) are clearly different. 
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Fig 

as the 



To ensure that the abundances are really different, 
the Kolmogorov-Smirnov test was used. This test gives 
the probability that two distributions are drawn from the 
same parent population. 



It was already shown in North et al. ( 1998 ) that the 
log g distribution of metallic and normal giants are differ- 
ent with a probability of 99.99%. However, the T c g distri- 
butions are comparable. The distributions of Tasini and 
Cmicro estimated by our method are also comparable. 

Fig. [ll] shows the cumulative distribution of abun- 
dances for each of the 6 elements that have the largest 
differences. They are Ca, Ti, Mn, Fe, Ni and Ba. Al and 
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Fig. 10. Abundance pattern of normal giants (■) and of metallic giants (A) 



Cr also have different distributions. Sc is not added, even 
though the corresponding distributions are clearly differ- 
ent, because of its correlation with V sini. 

Thus, we conclude that metallic giants have abun- 
dances of Al, Ca, Ti, Cr, Mn, Fe, Ni and Ba that are 
significantly larger than those of normal giants. 

7.3. Abundance variation with evolution 

Recently, the Montreal group has published a new study 
for normal stars. Models including all effects of atomic dif- 
fusion and radiative accelerations have been pres ented for 
stars of diffe rent masses in Turc otte e t al. ( 199£ ), Richer 
et al. fl2000| ) and Richard et al. ( [200 

The models of Turcotte et al. predict the evolution 
of the abundances for an F star of a given mass consis- 
tently with the internal structure. They are computed in 
the diffusion framework with only the mixing length as 
free-parameter. They include the effect of 28 chemical ele- 
ments (Z < 28) but no macroscopic mixing. Atomic diffu- 
sion plays an important role for stars more massive than 
1.3 M Q . We agree with Varenne & Moni er (1999 ) that the 
theoretical prediction of Turcotte et al. ( 1998 ) for F stars 
are not confirmed by observations. In our case, abundance 
of C, O, Si, Ca, Fe, and Ni do not show any variations with 
effective temperatu re (se e Fig. [l2| ). 

Richard et al. ( |2001 ) present stellar models that are 
evolved from the pre-main sequence to the giant branch for 
stars of 1.3 to 4.0 M Q and metallicities of Z=0.01 to 0.03. 
In order to compare our observations with these models, 
our sample was split into 7 mass groups with masses below 



1.4 M Q , 1.4 to 1.6 Mpi, 1.6 to 1.8 M P 



0; 



*0> 



2.4 to 2.6 M P 







and above 2.6 M Q . Fig. |13j shows the behaviour of Ca, Fe 
and Ni for stars of 1.6 to 1.8 M Q and 1.8 to 2.0 M Q with 
respect to logg. It can be directly compared with Fig. 10 
of Richard et al. (2001). They predict a large abundance 
variation of at least 1 dex with decreasing log g (or equiv- 
alently, with time ). In our case, there is no variation with 
\ogg except for stars of 1.8 to 2.0 Mq, and even in this 
case, the amplitude is much smaller than the theoretical 
predictions. In order to confirm the statistical significance 
of the correlation, 2x2 contingency tables have been com- 
puted, segregating stars with gravities larger or lower than 
\ogg = 4.0 and with [Fe/H] or [Ni/H] larger or lower than 
0.0. The corresponding % 2 were co mpute d taking Yates' 
correction into account (Crow et al. 196fj|) and the results 
were 8.98 and 10.94 for [Fe/H] and [Ni/H] respectively, 
which comfortably exceeds the 7.88 value corresponding 
to a 0.5% probability that abundances and logg are un- 
related. Thus, taken at face value, the correlation is sta- 
tistically significant. If this variation is physically true, it 
suggests that additional hydrodynamical processes are at 
play in these stars, which lower the effects of microscopic 
diffusion. However, before concluding to the physical real- 
ity of this correlation, it is necessary to test how far this 
behaviour can be explained by a correlation of abundances 
with Vsini or by NLTE effect. 



Takeda & Sadakane (|1997| ) also found a variation of 
iron abundance with logg, which was qualitatively the 
same, for Hyades' members. They explain this variation 
as a superficial phenomenon due to the existence of a mu- 
tual relationship between logg and V^sini. They found an 
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Fig. 11. Cumulative distributions of abundances of metallic giants (solid line) and of normal one (dashed line). The 
number in each panel indicates the probability that both distributions are drawn from the same parent distribution. 



anticorrelation between Vsini and logg, in the sense that 
high-gravity stars tend to be slow rotators. They consider 
that the cause of this relation is rotation, by which the 
apparent and effective values of log g are affected. In their 
case, log g is estimated from Stromgren photometry, and 
rotation affects this determination: it appears that logg 
values progressively decrease with increasing V smi. 

However, in our case, log g values were estimated from 
HiPPARCOS data, and they are not correlated with Vsinz 
(see Fig. |l4|). The effect of rotation is to increase the 
B2 — VI index. Maeder (1971) has shown that, to a good 



approximation, this index is essentially sensitive to the 
product Vs'mi. For a star of 1.4 M (2.0 M Q respectively), 
a Vsmi of 200 km s _1 leads to a decrease in log(T c ff) 
of 0.016 (0.010 respectively). The effect on M v is of a 
few hundredths of a magnitude (~ 0.1 mag respectively) 
(Maeder & Peytremann, 1970| ). This affects slightly the 
mass estimate when interpolating through evolutionary 
tracks of non-rotating models, but the error cannot ex- 
ceed 0.1 M Q (5% at 2 M ), propagating into no more 
than ~ 0.02 dex in log g. In fact our estimates of log g are 
close to the effective gravity, i.e. the sum of gravitational 



and centrifugal accelerations average over the visible disk; 
we verified that they agree quite well with photometric 
estimates computed from the uvbyfl colours and the cal- 



ibration by Moon & Dworetsky (1985). The interesting 
thing with effective logy is that for Vsmi ~ 200 km s _1 , 
it differs by no more, on the equator, than 0.1 dex from 
a non-rotating 2 M© star near the ZAMS, while the dif- 
ference is at least twice as large at the turn-off. Thus the 
observed logg axis (e.g. in Fig. 13) is slightly stretched 
for fast rotators, relative to the logg non-rotating stars 
would have in the same mass range. But since slow ro- 
tators are not affected, the net effect is just to blur very 
slightly the relation metallicity vs. logg shown in Fig. [l3| 
Therefore, this cannot explain the beha viour of \^r] with 
logg. Moreover, in Takeda & Sadakane (1997), logg varies 
on a range that is no more than 0.2 dex; under the assump- 
tion that this variation is due entirely to a bias in the log g 
estimate, it can hardly explain a 0.4 dex variation of [^1 
for the single line of neutral iron they used. 

Given the variation of T e g and logg for these stars, 
NLTE effect s might play an important role. Rentzsch- 
Holm (1996) has estimated NLTE effects on iron abun- 
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Fig. 12. Abundances of dwarf F-type stars as a function of T e g. 



dances for A-type stars. NLTE correction increase with 
r Gff go ing fr om 7000 to 10 000 K (see Fig. 4 of Rentzsch- 
Holm ( |1996D ). On the other hand, at given effective tem- 
perature, they increase with decreasing gravity. Since both 
T e ff and log g decrease with time when a star is evolving on 
the main sequence, this tends to flatten the NLTE correc- 
tion. For our stars of 1.8 to 2.0 M Q , there is a small NLTE 
correction of 0.1 to 0.15 dex. However, this correction is al- 
most constant (variation of 0.05 dex) for the T e g and logg 



considered^. Moreover, with our method, we do not have 
separate abundance values for Fe I and Fe n. Let us recall 
that Fell lines are almost not affected by NLTE effects. 
When T e ff increases, the ratio of line numbers of Fei and 
Fen decreases. Therefore our determination becomes less 
sensitive to NLTE effects. This also tends to flatten the 
variation of NLTE correction with T c ff . Therefore, we con- 



Note that this is valid only if T e ff and log g are correlated. 
The Pearson product-moment correlation coefficient is 0.972 
(0.968) for stars of 1.8 to 2.0 M (1.6 to 1.8 M respectively), 
and in the same order for the other mass interval. 
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Fig. 13. Relation between abundances of Ca, Fe and Ni with logg. Left: stars of 1.6 to 1.8 M . Right: stars of 1. 
to 2.0 M Q . The broken lines show the limits of the 2x2 contingency tables mentioned in the text. 



elude that, for stars of 1.8 to 2.0 M Q , only a small NLTE 
correction of about 0.1 dex, that is roughly constant with 
T e g, has to be taken into account when iron abundance 
is estimated with our method. To summarize, all points 
in the vs logg diagram for stars of 1.8 - 2.0 M Q in 
Fig. |l3| should be shifted upward by about 0.1 to obtain 
a more realistic picture. 



Thevenin & Idiart ( 1999 ) discuss the problem of log g 
determination when dealing with NLTE correction. Their 
work focuses on cooler stars, but the problem of logg de- 
termination is comparable. When the NLTE correction 



is applied to LTE abundances, abundances of Fei and 
Fen are no longer the same (assuming that logg had 
been adjusted to satisfy the LTE ionization equilibrium). 
Therefore, logg has to be corrected until ionization equi- 
librium is satisfied after NLTE correction. Nissen et al. 



( |1997| ) have shown that there exists a discrepancy between 
spectroscopic logg (adjusted assuming LTE) and those 
deduced from HiPPARCOS parallaxes. Thevenin & Idiart 
( 1999Q have solved this discrepancy by adjusting log g after 
NLTE correction; the latter are in close agreement with 
those deduced from HiPPARCOS (our logg are estimated 
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Fig. 7. Behaviour of Sr, Sc, Na, Fe with respect to Vs'mi. 



in this way). However, in order to apply NLTE corrections 
precisely, it is necessary to have separate abundance es- 
timates of Fei and Fen, which is not our case. Roughly, 
our determination can be seen as an average of both abun- 
dances weighted by the number of lines of each ionization 
stage. 

To summarize, NLTE effects on iron is no more than 
0.1 dex and remains constant to better than 0.1 dex, for 
stars of 1.8 to 2.0 M© when abundance are estimated with 
our method. 



Fig. 9. HR diagram of giant s tars w ith some some evolu- 
tionary tracks (Schaller et al., 1992 ). ■ represent normal 
giants with Am.2 < 0.013 and A represent metallic giants 
with Am 2 > 0.013. 



3.B 



Fig. 14. Plot showing V smi as a function of log g for stars 
of 1.8 to 2.0 M Q . 



In view of the above discussion, it seems that the cor- 
relation between [^1 and logg for 1.8-2.0 M is indeed 
physical, and due neither to systematic errors nor to se- 
lection effects. It agrees qualitatively with the prediction 
of Richard et al. (2001), but it is strange that only the 
stars of 1.8 to 2.0 Mq follow their models. It is true that 
more massive stars are less numerous, so their statistics 
is less good. For less massive stars, models predict a vari- 
ation that occurs on a shorter log g range. It is possible 
that this variation is too quick to be observed, but then 
we would expect to see at least an increase of the scatter 
for large logg values, while this is not the case. 

The comparison with the Montreal models shows that 
it is necessary to add hydrodynamical processes in or- 
der to reproduce observations, like mass loss and merid- 
ional circulation, as already suggested by the Montreal 
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group. Maybe turbulence will also have to be modelled in 
a slightly different way yet to be defined. 

We also note that the metal-rich 1.8 to 2.0 M Q stars 
lie just around the blue border of the instability strip or 
outside of it: the coolest of these (HD 144197) has T cS 
= 7691 K and the next cooler one (HD 148367) has T ofr 
= 7962 K, while the blue limit of the instability strip is 
defin ed by B2 - VI = -0.01 on the ZAMS (Jasniewicz 
1984) which correspond to T e g = 8016 K according to the 



calibration of Hauck & Kunzli (199£). The hottest metal- 
rich star (HD 195943) has T cff = 9057 K. Almost all stars 
with [Fe/H] < 0.0 are cooler and lie, therefore, either in- 
side the instability strip or to the red side of it. This might 
be more than a coincidence, and is reminiscent of the well 
known exclusion between S Scuti type pulsations and the 
Am phenomenon. The difference is that while most Am 
stars do lie inside the instability strip even though they do 
not pulsate, the metal-rich 1.8 to 2.0 M Q stars are found 
either outside or on the verge of it. 

8. Conclusion 

The method of analysis presented in Paper I was used to 
derive detailed abundances of 140 A and F-type stars with 
V smi up to about 200 km s . Vshii derived with this 
method are consistent with those in the literature. Only 
Cr, Sc, and Na are very sensitive to rotational velocity in 
the sense that their abundance determination is altered 
when Vsmi > 100 km s _1 , for various reasons. 

This study allows us to specify more exactly the na- 
ture of metallic A-F giants. They are more evolved than 
normal giants and they show overabundances of at least 
8 elements with respect to normal giants. These elements 
are Al, Ca, Ti, Cr, Mn, Fe, Ni and Ba. 

Our results do not confirm, in general, the predictions 
of the Montreal models. However, stars of 1.8 to 2.0 M Q 
show a clear correlation of iron abundances with logg 
which is reminiscent of theoretical predictions by Richard 
et al. (2001). Neither rotational velocity nor NLTE effects 
can explain this behaviour, which appears, therefore, to be 
real. The observed [4^1 and \^jf] variations qualitatively 
agree with those predicted by the Montreal models, but 
are much weaker. The smaller amplitude indicates that 
their models need more hydrodynamical processes (merid- 
ional circulation, mass loss, . . . ) and/or a better tuning 
of convection. Our large homogeneous set of abundances 
hopefully can help this work. 
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